Hole-pair hopping in arrangements of hole-rich/hole-poor domains in a quantum 

antiferromagnet 



(N 
O 

o 

> 
O 

in 

(N 



c3 



O 
o 



> 

O 

in 



(N 

o 

c3 



S3 
O 

o 



X 
S3 



Efstratios Manousakis 

Department of Physics and MARTECH, Florida State University, Tallahassee, Florida, 32306, USA and 
Department of Physics, Solid State Physics Section, University of Athens, Athens, Greece 

(Dated: February 1, 2008) 

We study the motion of holes in a doped quantum antiferromagnet in the presence of arrangements 
of hole-rich and hole-poor domains such as the stripe-phase in high-Tp cuprates. When these 
structures form, it becomes energetically favorable for single holes, pairs of holes or small bound-hole 
clusters to hop from one hole-rich domain to another due to quantum fluctuations. However, we find 
that at temperature of approximately 1QQK, the probability for bound hole-pair exchange between 
neighboring hole-rich regions in the stripe phase, is one or two orders of magnitude larger than 
single-hole or multi-hole droplet exchange. As a result holes in a given hole-rich domain penetrate 
further into the antiferromagnetically aligned domains when they do it in pairs. At temperature 
T ~ lOOif and below bound pairs of holes hop from one hole-rich domain to another with high 
probability. Therefore our main finding is that the presence of the antiferromagnetic hole-poor 
domains act as a filter which selects, from the hole-rich domains (where holes form a self-bound 
liquid), hole pairs which can be exchanged throughout the system. This fluid of bound hole pairs 
can undergo a superfluid phase ordering at the above mentioned temperature scale. 

PACS numbers: 74.20.-z,71.10.-w,71.45.Lr 



I. INTRODUCTION 

When an isotropic quantum antiferromagnet on an in- 
finite square lattice, as described by the t — J model with- 
out the inclusion of the long-range part of the Coulomb 
interaction, is doped with two holes a bound state is 
formed B. This bound state is not a sign of a pairing in- 
stability but that of a phase separation instability B |j or 
possibly stripe ordering Q. While in the high-T c cuprates 
stripes have been experimentally observed, at the theo- 
retical level, there are two different scenarios to explain 
their formation. In the first scenario, when the antifer- 
romagnet is doped with more than two holes they form 
a larger droplet and there is evidence || ||] that even for 
relatively small values of J/t there is a critical electron 
density n ps (J/t) above which the system is separated into 
two phases: one phase in which the holes are bound into 
a hole-rich region and an all electron phase characterized 
by antiferromagnetic order. When one adds the long- 
range part of the Coulomb interaction this phase is for- 
bidden because of the fact that the energy density of a 
droplet with finite charge density blows up as the droplet 
size goes to infinity. The system can accommodate this 
tendency for phase separation by forming a state which 
is made out of alternating microscopic regions of the two 
phases but with a zero net charge when averaged over a 
larger region. Such possible states are shown in Fig. [ljand 
in Fig. |[ In the second scenario the stripe-ordered state 
is obtained right away as the ground state of the t — J 
model without the inclusion of the long-range Coulomb 
interaction j|. The validity of these calculations has been 
discussed^ and the issue has not been resolved unequiv- 
ocally yet. 

In this paper we consider the electronic system in a 
state characterized by alternating regions or domains 



with different charge such as the static or dynamic stripe 
phase (Fig. [I]) or that of Fig. || or other structures which 
take care of the tendency for phase separation and the ef- 
fects of long-range Coulomb repulsion; when such states 
form, it becomes energetically favorable for single holes 
or small bound-hole clusters to hop from one hole-rich 
domain to another due to quantum fluctuations. We 
show that the probability of bound hole-pair hop is one 
to two orders of magnitude larger than that of single-hole 
or multi-hole droplet hop between neighboring hole-rich 
regions. Therefore, pairs of holes in a given hole-rich do- 
main prefer to penetrate further into the antiferromag- 
netically ordered domains relative to single or multi-hole 
clusters and they have a much higher chance to tunnel 
from one hole-rich domain to another. In particular, we 
find that below T ~ 100 A' these bound hole-pairs tun- 
nel through the antiferromagnetic domains at high rates. 
Therefore, the system can be pictured as consisted of two 
inter-penetrating subsystems a) a subsystem of hole-rich 
domains in which the holes exist as a self-bound liquid 
and b) a fluid of bound hole pairs which exists in the anti- 
ferromagnetically ordered hole-poor domains. The latter 
fluid is composed of bound hole-pairs with a large char- 
acteristic binding energy scale of the order of the phase 
separation energy or the stripe formation energy scale. 
The main point of this paper is that the existence of the 
antiferromagnetic domains act as a filter which allows 
only pairs to exist over the entire system. 

The second part of this paper includes a discussion 
which involves some speculation. The case of strong pair- 
ing correlations within each hole-rich strip and the role 
of the inter-strip coupling has been examined by other 
authors^). We wish to consider the limit where pairing 
within each hole-rich domain is weak at temperature of 
the order of 1Q0K. We argue that the fluid of bound 
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hole-pairs, which is selected by the presence of the inter- 
vening antiferromagnetic regions and exists throughout 
the system, can undergo a superfluid phase ordering at a 
relatively high critical temperature. In the limit of weak 
pairing within the hole- rich domains, this critical temper- 
ature associated with superfluid phase ordering depends 
on the pair effective mass inside the antiferromagnetic 
domain and the distance between the hole- rich domains. 



II. SINGLE-HOLE, HOLE-PAIR AND 
MULTI-HOLE TUNNELING 



In order to understand single-hole, hole-pair or multi- 
hole tunneling from one hole-rich domain (which could 
be a hole-rich strip or droplet) we will use a continuum 
model of holes of effective mass m* , interacting through 
a mean field V e ff introduced by the environment and 
the stripe-ordered or the dynamic-stripe configuration 
and via a residual hole-hole interaction V. The partition 
function for the system of N holes can be written, using 
Feynmanpl no) path integral formulation, as follows: 
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where the sum is over all N\ hole-permutations P and 
[P] is the order of permutation. The integration is over 
paths Xi(r) in imaginary time r. The final configuration 
x*i((i) at r = (3 is any permutation Pxi(0) of the initial 
hole-positions. 

Let us consider a state of alternating hole-rich and 
antiferromagnetically aligned hole-poor domains. When 
these alternating regions order they form the stripe- 
ordered state shown in Fig. [l](a). In addition we will 
keep in mind other inhomogeneous configurations such 
as the one in Fig. || In the case of the stripe state, each 
hole experiences a potential V e f / with a profile shown in 
Fig. |l](b). We wish to understand the tunneling of one 
or more holes from one hole-rich domain to the nearest 
hole-rich domain. First, let us consider the motion of a 
single- hole and estimate the factor exp(—S e f /) for a hole 
displacement from the boundary of a hole-rich domain by 
a distance x within imaginary time r inside the antifer- 
romagnetic hole-poor domain. At the beginning the hole 
is near the boundary of the hole-rich with the hole-poor 
region and so an estimate of the S e f t is given as follows 



S$ f (x,T) 



2h 2 t 



Vir. 



(2) 



Here ml is the hole effective mass inside the antiferro- 
magnet (hole-poor region), V\ is the difference between 



the hole energy when the hole is in the hole-poor region 
and the hole energy when the hole is in the hole-rich re- 
gion. This is determined by the presence of an effective 
barrier imposed by the existence of the two different do- 
mains and more significantly by the hole-hole interaction. 
A single hole sticks to the hole-rich region strongly while 
a pair of holes sticks to the rest of the hole-rich region 
relatively weakly. This is shown on the basis of the t — J 
model later below. This means that the value of the bar- 
rier V\ for a single hole to enter the antiferromagnetic 
region is significantly larger than the barrier Va which is 
experienced by a pair of holes. The path with the most 
significant contribution to the partition function is that 
which corresponds to a characteristic time scale r (the 
time for which the hole stays in the hole-poor region) 
such that S^(x, r) for a given x is minimum. The min- 
imum of the action (M) is obtained for 
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and for this value of r the action is given by 
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This gives the well-known result for the quantum me- 
chanical penetration probability (e _s< ). Therefore the 
"penetration" depth inside the hole-poor region (the 

value of x where e _s< ' = 1/e) is given by 



Ai = 



(5) 



As we will show below a pair of holes sticks much less 
strongly to the hole-rich region; thus, the difference in 
energy V2 for the pair in the hole-rich and hole-poor do- 
mains is much smaller than V\. The penetration depth 
for hole-pair or multi-hole tunneling is given by 
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We would like to estimate these penetration depths and 
probabilities on the basis of results obtained by Green's 
function Monte Carlo simulation of the t — J model jlj ||, 
0,0- 

We consider the t — J model in the parameter range 
where phase separation was found H, [|, ^2). Starting 
from the undoped insulator we introduce N holes in a 
system of N s — L x L sites. The ground state of the 
system keeping N finite and L — > 00 is expected to be 
a two-dimensional droplet with energy per hole given by 
e(N) = (E(N) - E(0))/N. The situation is shown in 
Fig. I for J/t = 1 and in Fig. | for J/t = 0.4. The circles 
with the error bars give the energy per hole e(N) as a 
function of the number of holes. The solid line in both 
figures denotes the minimum of the energy per hole as 
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FIG. 1: (a) A striped ordered state, i.e., alternating strips of 
the hole-rich and the hole-poor phases. Two successive hole- 
poor domains are separated by a pi-phase shift. In general we 
consider the case of dynamic stripes where while these alter- 
nating domains exist there is no stripe ordering, (b) The ef- 
fective potentials experienced by a single hole (solid line) and 
bound hole-pair when attempting to tunnel from one hole- rich 
region to another. The corresponding wave functions of the 
single-hole and bound hole-pair are also shown schematically 
as ipi and ip2- 



0.5 



0.4 



° 0.3 - 



0.2 



0.1 - 



J/t = 1 



Minimum 



1 2 3 4 5 6 7 
Number of holes 



FIG. 3: The energy per hole of a hole-rich droplet in an an- 
tiferromagnetic background as a function of the number of 
holes at J/t — 1 (dashed line and solid circles). The solid line 
is the value of the minimum value of the energy per hole as a 
function of the hole density for J/t = 1 
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FIG. 2: A hexagonal superlattice formed from droplets of the 
hole-rich phase inside an antiferromagnetic background. 



a function of the hole density, namely the energy of the 
phase separated state e ps . The energy per hole e(N) is 
always higher than e ps . The quantity Vn discussed pre- 
viously, is given by (e(N) — e ps )N. The important facts 
to notice is that a) the energy per hole for the bound-pair 
of holes is quite close to that of the energy per hole in the 
phase separated state and b) introducing two holes inside 
the quantum antiferromagnet lowers the energy per hole 
significantly relative to the single hole case and further 
introduction of holes leads to formation of droplets with- 
out as significant gain in binding energy. This implies 
that it is far easier to separate a hole-pair than a single- 
hole from the hole-rich region. The energy of a hole-rich 
droplet inside a quantum antiferromagnet has an addi- 
tional contribution due to long-range Coulomb repulsion. 
It is clear than when we switch on this interaction, the 



FIG. 4: The energy per hole of a hole-rich droplet in an an- 
tiferromagnetic background as a function of the number of 
holes at J/t = 0.4 (dashed line and solid circles). The solid 
line is the value of the minimum value of the energy per hole 
as a function of the hole density for J/t = 0.4 



macroscopically phase separated state cannot be realized 
because the energy per hole has positive infinite energy. 
However, a state of microscopic phase separation is al- 
lowed because there is no net charge on a large scale and 
examples of such states are given in Fig. [j] and Fig. ||. 
The effect of this long-range interaction does not alter 
our conclusion and it will be discussed later. 

While m% increases almost linearly with the number 
of holes ArjlJ, the function V N = (e(N) - e ps )N drops 
sharply at the value N = 2. As a result the combina- 
tion of m* N VN which occurs in Eq. ^ causes a peak in Ajy 
at N — 2. The fact that the largest Xn corresponds to 
bound hole-pairs implies that when stripes form, because 
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of the overlap of the wave- function of hole clusters within 
each hole-rich region with the wave-function of a neigh- 
boring hole-rich strip, the number of bound hole-pairs 
which will be shared among all strips is much larger than 
the number of single holes or the number of any other 
multi-hole clusters shared by all the strips. In addition 
using values for Vjv and m%/m obtained from numerical 
studies of the t — J model 0, |, |ll|. ||] and the Hubbard 
model|L| we find that A2 is of the order of 10A which 
is of the inter-strip distance when the stripes order. In 
the next section we show that the contribution of such 
tunneling configurations (in which hole-pairs cross from 
a hole-rich strip to a neighboring hole-rich strip through 
the intervening antifcrromagnctic hole-poor domains) to 
the path integral (Q) become very significant at and below 
a temperature of approximately 100A. 

Tunneling of hole-pairs between neighboring strips 
contribute to the path-integral (Q) by a factor 
exp{—S pa i r {d, (3)) while when a single-hole crosses from 
one hole-rich strip to a neighboring one through the anti- 
ferromagnetic domain it contributes to the path integral 
by an amount exp{—S\{d,j3)). The ratio of these contri- 
butions is given by exp(— AS , i 2 (d, 0)) where A5i2(d, (3) 
can be estimated as 

AS 12 (d,p) = ^=^j + (V! V 2 )(3. (7) 

Taking T = 100 A and using the available results for 
J/t = 0.4, namely V x ~ 0.07i, V 2 ~ 0.02£, and the 
values of rrit / mp ~ 6.5,7712/^0 ~ 22.5 (as obtained by 
Trugman(l3) for U/t = 10) and the value d ~ 16A (a 
typical distance between strips in a stripe-ordered state) 
we find that Si — S pa i r ~ 4. This implies that the contri- 
bution of a single hole e~ Sl is approximately 1-2 orders 
of magnitude smaller than that of a bound hole-pair at 
the above mentioned temperature. 

In order to give an accurate determination of such tran- 
sition rates for hole-pairs we need a) a more accurate 
determination of the pair binding energy difference rela- 
tive to the energy of the phase separated state and b) to 
know the effect of the long-range Coulomb repulsion on 
the binding energy of two holes inside an antiferromag- 
netic domain relative to the energy of these two holes 
when they are in the hole-rich domain. The value we used 
above, namely V 2 ~ 0.02t will be significantly reduced in 
the presence of the long-range (LR) Coulomb interaction. 
The turning on of the LR Coulomb interaction will reduce 
the binding energy of a bound hole-pair to the rest of the 
hole-rich region. The stripes result from the competition 
between the tendency of the system for phase separation 
and the LR Coulomb interaction. Therefore the hole-rich 
region contains just enough holes to keep the balance be- 
tween these two opposing tendencies. This implies that 
the binding energy of a bound pair to the rest of the 
hole-rich domain relative to the case where the hole-pair 
is inside the antiferromagnetic hole-poor region is very 
small. This is in agreement with the findings of Arrigoni 
et al. [M who studied the t — J model and they added the 



LR Coulomb interaction using a combination of density 
matrix renormalization group method for the short-range 
part and a Hartree approximation to take into account 
the LR part. When they added the LR Coulomb inter- 
action they find enhanced superconducting correlations 
which was associated to the enhancement of pair tunnel- 
ing between stripes. 

III. SUPERFLUIDITY OF THE GAS OF 
BOUND HOLE-PAIRS 

This section includes a discussion which is speculative. 
If pairing correlations within a given hole-rich stripe do- 
main are strong we are in a limit which has been consid- 
ered by other authors Q. In such case a weak inter-strip 
coupling can give rise to a high Tq and superconductiv- 
ity competes with charge density wave ordering. Here 
we wish to examine the case where pairing correlations 
within a given strip are weak at temperature of the or- 
der of 100K. Strong pairing correlations exist in ladders 
but when stripes occur in a two-dimensional system it is 
uncertain that such one-dimensional pairing correlations 
are strong at this temperature. 

If we neglect V2, S pa i r ~ (m^d 2 /2H 2 )KbT and it is 
of order of unity at T ~ 100 A . At around this temper- 
ature the system of such hole-pairs acquires long-range 
phase coherence and non-zero winding number due to 
long chains of exchanges of hole-pairs between neighbor- 
ing strips. The most significant contributions to the par- 
tition function (Q) at temperature of the order of 100A 
and below can be written as 

z = ^ (2^9) 2 /n^.a-^) 

where N p = N/2 is the number of hole pairs in the 
system with a distribution governed by the function 
g(zi,Z2,...,ZN p ) giving the probability for a given con- 
figuration {zi} of the bound hole pairs as in a Bose 
fluid. There is a sum over all possible permutations 
of hole-pairs. Here we have assumed that the path in- 
tegration over all trajectories in Euclidean-time in the 
path-integral of Eq. (Q) has been carried out(T(J]. Let 
us consider a typical displacement of a hole inside the 
hole-rich domains which is associated with the factor 
fi = cxjp(—m' 1 \x\ 2 /2h 2 /T — U\{t)) where U\(t) is the 
interaction part of the action and m\ the hole effective 
mass inside the hole-rich domains. Similarly when a 
pair of holes moves a similar typical distance between 
pairs inside the hole-rich domain we can associate a fac- 
tor f 2 = exp{~m 2 \x\ 2 /2h 2 /T - U 2 (t)). For hole dis- 
placements of the order of the inter-hole distance in the 
hole-rich domains, both fi and f 2 become of order of 
unity for significantly smaller values of the imaginary 



5 



time r (higher temperatures) as compared to the time 
scale required for hole-pair exchange between neighbor- 
ing hole-rich domains. Namely, in order to provide an 
optimum hole distribution within the strips the holes, 
in these regions, adjust their positions in shorter times 
scales; the effects discussed previously which involve tun- 
neling of pairs of holes between two neighboring hole- 
rich regions are longer-time effects. Therefore we have 
assumed that at such lower temperature where these lat- 
ter effects become energetically favorable, we can inte- 
grate out these short times scale configurations; the holes 
within each strip have enough time to relax to a distri- 
bution g(z\, ...,Zk) which is experienced by the bound 
hole-pairs as they move from one hole-rich domain to an- 
other. This distribution g describes the stripe ordering 
and the charge-stripe correlations. We have eliminated 
the single-hole degrees of freedom even within a strip. 
We consider the hole-rich domains as made of pairs of 
holes with the same charge distribution because what is 
important in describing our ideas is the movement of the 
gas of bound hole pairs in the space between hole-rich do- 
mains. The space occupied by the hole-rich strips plays 
only the role of a reservoir from where these bound hole 
pairs emerge. The characteristic time-scale associated 
with the motion of these degrees of freedom in imagi- 
nary time is much shorter than the inverse temperature 
at which tunneling of bound hole-pairs between hole-rich 
domains becomes likely. 

When a bound hole-pair is exchanged between 
neighboring stripes, it leads to a factor / = 
exp(— ■m\j(2fih 1 }{zi — zj) 2 ). This can occur anywhere in 
the system between any two neighboring stripes. These 
pair-hole movements cannot occur in an uncorrelated 
manner because when the pairs move to another hole- 
rich domain they disturb the balance of charge. Such 
unbalanced configurations correspond to large values of 
the effective action and they will not contribute to the ob- 
servables. When closed rings of exchanges occur, namely 
cyclic permutations of bound-pairs, any transient charge 
imbalance should be eliminated. Below a certain temper- 
ature, of the order of 1QQK as shown above, the factor / 
becomes of the order of unity and then such coherent mul- 
tiple exchanges can occur which lead to the off-diagonal- 
long-range order. The arguments which lead to this last 
part of the picture are very similar to those discussed by 
FcynmanjlO for the case of liquid helium-4. 

Since the partition-function (||) describes a system of 
purely boson degrees of freedom, we can gain additional 
support for this picture by discussing what happens in 
a simpler model. We consider a system of bosons in- 
teracting via a Van der Waals interaction and we add 
a long-range a/r Coulomb- like interaction with a very 
small value for a. First let us consider the case of 
a = 0, namely that of a pure self-bound boson liquid 
on a two-dimensional substrate. Such a system has a 
physical realization, namely, it describes liquid helium-4 
on a substrate such as graphite. The pure liquid neu- 
tral helium-4 on a 2D smooth substrate (without sub- 



strate corrugations) with periodic boundary conditions 
has been studied Q using path- integral Monte Carlo 
simulation pl|. It was found that for values of the 2D 
density below the 2D equilibrium density and at low 
temperature the system forms a 2D liquid droplet with 
non-zero winding number (superfluid density) for densi- 
ties well inside the phase separated region. Let us discuss 
what we expect to find if we introduce the long-range a/r 
interaction term with a small value of a so that a self- 
bound droplet with at least two atoms exists as a result of 
the competition of the long-range and the Van der Waals 
attraction. We expect that the uniformity of the system 
will be spoiled by the presence of this long-range inter- 
action and droplets of certain size will be formed. The 
system optimizes the number of atoms Nq in a droplet, 
namely, the binding energy per atom will be very small 
compared to the magnitude of the van der Waals inter- 
action. This implies that an atom in any charged droplet 
will be only very weakly bound to the droplet and this 
will allow a non-zero winding number caused from ex- 
changes of atoms between these droplets. An estimate of 
the temperature scale where such exchanges occur is 



K B T C 



2md 2 



(9) 



where d is the average distance between neighboring 
droplets. We expect to find a similar result if one in- 
troduces a weak external potential so that the charged- 
superfluid is confined in a charge-stripe-like configura- 
tion; namely, the off-diagonal-long-range order sets in at 
a temperature of the above mentioned magnitude where 
d is the distance between successive hole-rich domains. 



IV. CONCLUDING REMARKS 

We have considered the state in which alternating hole- 
rich and hole-poor (antiferromagnetically ordered) do- 
mains form in a strongly correlated electron system. 

Here we have used existing numerical results on the 
t — J model to show that the presence of antiferromag- 
netically ordered regions between such hole-rich domains 
acts as filter which allows only pairs of holes to jump 
from one hole-rich domain to a nearest neighboring one 
at high rates. This implies that effectively there is a fluid 
of bound hole pairs which can move through the entire 
system. 

We further discussed a speculative point of view. We 
considered the limit in which pairing correlations within 
each hole-rich domain are weak which implies that such 
a domain is characterized by a small superconductivity 
ordering temperature scale. In this case we argued that 
because of the existence of the previously discussed fluid 
of bound hole-pairs, the characteristic energy scale as- 
sociated with phase coherence between such pairs is de- 
termined by the distance between neighboring hole-rich 
domains and the pair effective mass inside the antiferro- 
magnetically ordered domains. 
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